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ABSTRACT 



> . 

Aims. To investigate the possible relationships between dynamical status and other important characteristics of galaxy 
clusters, we conducted a study of X-ray cluster morphology using a sample of 101 clusters at redshift z~0.05-l taken 
from the Chandra archive. 

£ — ■ Methods. The X-ray morphology is quantitatively characterized by a series of objectively measured simple statistics of 

the X-ray surface brightness distribution, which are designed to be robust against variations of image quality caused 
by various exposure times and various cluster redshifts. Using these measures, we quantitatively investigated the rela- 
tionships between the cluster X-ray morphology and various other cluster characteristics. 

Results. We found: (1) Our measures are robust against various image quality effects introduced by exposure time 
difference, and various cluster redshifts. (2) The distorted and non-distorted clusters occupy well-defined loci in the 
L-T plane, demonstrating the measurements of the global luminosity and temperature for distorted clusters should be 
interpreted with caution, or alternatively, a rigorous morphological characterization is necessary when we use a sample 
of clusters with heterogeneous morphological characteristics to investigate the L-T or other related scaling relations. 
(3) Ellipticity and Off-center show no evolutionary effects between high and low redshift cluster subsets, while there 
may be a hint of weak evolutions for the Concentration and Asymmetry, in such a way that high-z clusters show more 
distorted morphology. (4) No correlation is found between X-ray morphology and X-ray luminosity or X-ray morphology 
and X-ray temperature of clusters, implying that interaction of clusters may not enhance or decrease the luminosity 
or temperature of clusters for extended period of time. (5) Clusters are scattered and occupying various places in the 
SpH, plane composed of two X-ray morphological measures, showing a wide variety of characteristics. (6) Relatively strong 

q . correlations in Asymmetry-Concentration and Offcenter-Concentration plots indicate that low concentration clusters 

generally show high degree of asymmetry or skewness, illustrating the fact that there are not many highly-extended 
smooth symmetric clusters. Similarly, a correlation between Asymmetry and Ellipticity may imply that there are not 
. many highly-elongated but otherwise smooth symmetric clusters. 
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■ 1. Introduction ies, and lensing. However, the distribution study requires a 

large number of galaxies, and is more susceptible to contam- 
Over the past decade, studies have provided evidence that a ination from foregr ound and background objects. Lensing is 
significant fraction of galaxy clusters have undergone recent alg0 sensitive to this contamination, and does not have good 
mergers (e.g. Geller & Beers 1982; Dressier & Shectman gpatial resolution except for the central region of a cluster . 



These mergers are observed as disturbed cluster An alternative method comes from x _ ray wave lengths, be- 
morphologies. The important connection between the mor- cauge dugter mergers compress and heat the intracluster 
phologies of galaxy clusters and cosmological parameters gaS; and this can be measured as distortions of the spa- 
has received much recent attention (Richstone, Loeb, & tial distribution of X -ray surface brightness and temper- 
Turner 1992; Evrard et al. 1993; Mohr et al. 1995). This ature Moreover, X-ray emissivity is proportional to the 
connection has generally been formulated m terms of the square of thc electron density, and therefore less affected 
frequency of 'substructure in clusters and from qualitative by the superposed struc tures than optical data. Jones & 
measures of the frequency of substructure m clusters, mves- Forman ( 1999 ) visua n y examined 208 clusters observed 
tigators have attempted to determine fi ro (e.g. Richstone with Einstein X -ray satellite and separated these clusters 
et al. 1992) and the power spectrum of primordial density intQ gix morphological classes . They found that about 40% 
fluctuations (e.g. David et al. 1993) by comparison to Press of their clusterg displayed some type of 'substructure'. 
& Schechter (1974) type predictions of the distribution of 

collapsed objects. However, a more quantitative measure of cluster struc- 

Methods to quantify structures at optical wavelengths ture at X-ray wavelengths is desirable to quantitatively 

have mostly used both the distribution of cluster galax- test various scenarios related to clusters, including cosmol- 

ogy Using Einstein images, Mohr et al. (1995) measured 

Send offprint requests to: Y. Hashimoto, emission-weighted centroid variation, axial ratio, orienta- 

e-mail: hashimot@mpe.mpg.de tion, and radial falloff for a sample of 65 clusters, while sev- 
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eral other studies used ellipticity (e.g. McMillan et al. 1989; 
Gomez et al. 1997; Kolokotronis et al. 2001; Melott et al. 
2001; Plionis 2002). Buote & Tsai (1995, 1996) used a power 
ratio method for 59 clusters observed with ROSAT, while 
Schucckcr et al. (2001) conducted a study of 470 clusters 
from ROSAT-ESO Flux-Limited X-ray (REFLEX) cluster 
survey (Bohringcr et al. 2000), using sophisticated statis- 
tics, such as Fourier elongation test, Lee test, and (3 test. 

Despite the success of these studies, all of them are un- 
fortunately limited to clusters in the nearby universe (z < 
0.3), where we may expect to see less frequent morphologi- 
cal distortions, and little evolutionary effect. This is due 
to the fact that, until recently, only a small number of 
high-z clusters have been known, or observed with sufficient 
depth and sufficient spatial resolution. With the advent of 
big-aperture satellites equipped with high spatial-resolution 
instruments, such as Chandra and X M M -Newton, to- 
gether with newly-available lists of distant clusters gen- 
erated based on various deep cluster surveys, it is finally 
possible to extend the morphological study to higher red- 
shifts. Indeed, Jeltema et al. (2005) have recently extended 
the power ratio method of Buote & Tsai to 40 clusters at 
z=0. 15-0.9 using Chandra data, and reported the evolution 
of cluster morphology between two rcdshift bins (z<0.5 and 
z>0.5). 

Extending the morphological study to high redshift is 
important but a difficult task because of inevitable low data 
quality associated with high-z clusters. Conventional meth- 
ods for characterizing the cluster X-ray morphology are of- 
ten sophisticated and some methods have an advantage of 
being more directly related to a particular characteristic 
of a cluster, such as mass, dark matter content, or gravita- 
tional potential. However, most of these conventional meth- 
ods are originally developed to analyze the low redshift clus- 
ters with high data quality, and, perhaps because of their 
intrinsic sophistication, ofter require many photon counts, 
making the measures rapidly uncertain or unmeasurable as 
the data quality decreases. Moreover, these methods also 
often require some interactive processes, and thus, are not 
suitable for the investigation involving a large dataset with 
a wide variety of image quality where various systemat- 
ics should be treated and removed in a consistent manner. 
Although it is important to try to extend the sophisticated 
methods to high redshift, a complementary study using the 
robust measures of the cluster morphology, less sensitive 
to variation in the data quality and suitable for a large 
dataset, is much needed, to enable us to study the low-z 
and high-z universe in a uniform manner. 

Here we report our study of X-ray cluster morphology 
using a sample of 101 clusters at redshift z~0.05-l taken 
from the Chandra archive. The X-ray morphology is quan- 
titatively characterized by a series of objectively-measured 
simple statistics of X-ray surface brightness, which are de- 
signed to be robust against variations of image quality 
caused by various exposure times and various cluster red- 
shifts. Using these measures, we quantitatively investigated 
the relationships between the cluster X-ray morphology and 
various other cluster characteristics. 

This paper is organized as follows. In sec 2, we describe 
our sample and data preparation, while in sec 3, details of 
our measures are described, and in sec 4, uncertainty and 
systematics are investigated. Sec. 5 summarizes our results. 
Throughout the paper, we use H Q — 70 km s" 1 Mpc -1 , 
f2 m =0.3, and SIa=0.7, unless otherwise stated. 



2. Sample & data preparation 

Almost all clusters are selected from flux-limited X-ray sur- 
veys, and data are taken from the Chandra ACIS archive. 
A lower limit of z = 0.05 or 0.1 is placed on the redshift to 
ensure that a cluster is observed with sufficient ficld-of-view 
with ACIS-I or ACIS-S, respectively. The majority of our 
sample comes from the ROSAT Brightest Cluster Sample 
(BCS; Ebeling et al 1998), and the Extended ROSAT 
Brightest Cluster Sample (EBCS; Ebeling et al. 2000). The 
BCS sample includes 201 clusters, with the flux limit of 
4.4xl0 -12 erg s _1 cm~ 2 (0.1-2.4 kcV). The authors esti- 
mated a sample completeness of 90 % for the 201 BCS 
clusters, and 75 % for the EBCS clusters. When combined 
with EBCS, the BCS clusters represent one of the largest 
and most complete X-ray selected cluster samples, and they 
are currently the most frequently observed by Chandra. 
As of 2005 October, 55 BCS + 13 EBCS (hereafter BCS) 
clusters with z > 0.05 (ACIS-I), or z > 0.1 (ACIS-S), are 
publicly available in the Chandra archive. Additionally we 
included all clusters from the X-ray flux limited sample 
of Edge et al. (1990) at z > 0.05 or 0.1 not in the BCS 
that were observed with the Chandra ACIS. This added 
12 more clusters. The Edge et al. sample is estimated to 
be ~90% complete, and contains the 55 brightest clusters 
from EX OS AT, HE AO - 1, and Einstein. 

To extend our sample to higher redshifts, additional 
high-z clusters are selected from various deep surveys; 10 of 
these clusters are selected from the ROSAT Deep Cluster 
Survey (RDCS: Rosati et al. 1998), 10 from the Einstein 
Extended Medium Sensitivity Survey (EMSS; Gioia et al. 
1990; Henry et al. 1992), 14 from the 160 Square Degrees 
ROSAT Survey (Vikhlinin et al. 1998), 2 from the Wide 
Angle ROSAT Pointed Survey (WARPS; Perlman et al. 
2002), and 1 from the North Ecliptic Pole survey (NEP; 
Gioia et al 1999), RXJ1054 was discovered by Hasinger et 
al. (1998), RXJ1347 was discovered in the ROSAT All Sky 
Survey (Schindler et al. 1995), and 3C295 has been mapped 
with Einstein (Henry & Henriksen 1986). 

The resulting sample we processed contains 120 clus- 
ters. At the final stage of our data processing, to employ 
our full analysis, we further applied a selection based on 
the total counts of cluster emission, (for details, please see 
Sec. 4), eliminating clusters with very low signal-to- noise 
ratio. Clusters whose center is too close to the edge of the 
CCD are also removed. The resulting final sample contains 
101 clusters with redshifts between 0.05 - 1.26 (median z 
= 0.226), and luminosity between 1.0 x 10 44 - 1.2 x 10 46 
erg/s (median 8.56 x 10 44 erg/s) (Fig.l). The final cluster 
sample together with their published redshifts and bolomet- 
ric luminosities, if available, as well as /3s and core radii, are 
listed in Table 1. 

We reprocessed the level=l event file retried from the 
archive using CIAO v3.1. and CALDB v2.29. For obser- 
vations taken in the VFAINT mode, we ran the script 
acis_process_events to flag probable background events, us- 
ing the information in a 5 x 5 event island. We also applied 
the charge transfer inefficiency (CTI) correction and the 
time dependent gain correction for ACIS-I data, when the 
temperature of the focal plane at the time of the observa- 
tion was 153 K. The data were filtered to include only the 
standard event grades 0,2,3,4,6 and status 0, then multiple 
pointings were merged, if any. We eliminated time intervals 
of high background count rate by performing a 3 a clipping 
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References: a: Reiprich & B6hringcr(2002); b: Ota & Mitsuda(2004); 
c: Ettori et al(2004); d: Kempner et al(2003); c: O'Hara ct al(2004); f: 
Ota et al. (2004) (R c & (3 from Ota & Mitsuda(2004); g: Hashimoto 
ct al(2004); h: David & Kempner (2004); i: Mazzota et al. (2001); j: 
Wu et al. (1999); k: Allen (2000); 1: Wang et al. (2004). Comments: 
(1): Distant southern component A115S is excluded; (2): Distant 
southern component A1758S is excluded; (3): RXJ0152, both north 
and south components are treaded as one cluster; (4): RXJ1054, 
both east and west components are treaded as one cluster. 



of the background level using the script analyzeJtcrv. To 
prepare the images for analysis, we selected photons in the 
observed-frame 0.7-8.0 keV and rest-frame 0.7-8 keV bands 
initially binned into 0.5 " pixel (see Sec 4 for the binning 
scheme of later analysis steps). We corrected the images 
for exposure variations across the field of view, detector 
response and telescope vignetting. 

We detected point sources using the CIAO routine 
cclldetect with signal-to-noise threshold for source detec- 
tion of three. An elliptical background annulus region was 
defined around each source such that its outer major and 
minor axes became three times of the source region. We 
removed the detected sources, except for a source at the 
center of the cluster which was mostly the peak of the sur- 
face brightness distribution rather than a real point source, 
and filled the source regions using the CIAO tool dmfilth. 
The images were then smoothed with Gaussian <r=5" . We 
have decided to perform the smoothing, as well as the total- 
count cut (see sec 4.2.3 for detail) to avoid the case where 
we have an image predominantly with zero count pixels, 
which makes the exposure map correction difficult, as well 
as the determination of the object region (see below), and 
for the investigation of various systematics (see sec 4). We 
found that the choice of smoothing-sigma hardly affects our 
robust morphological values (please see section 4.2.2. for 
detail) . 

Some clusters have a chip gap, or bad column inside the 
extracted cluster region. Most of these clusters, however, 
were observed with multiple pointings, thus those artifacts 
were reasonably corrected by exposure map. For those clus- 
ters with a single pointing, the artifacts were all crossing 
the cluster region far (typically more than 2 arcmin) from 
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Fig. 1. Distributions of the X-ray bolometric luminosity 
(Lbol) and redshift (z) for our sample. 



the cluster center. Using clusters with multiple pointing ob- 
servations, we discovered that, the effect of these artifacts 
are negligible on our morphological measures, particularly 
after the smoothing. 

We decided to use isophotal contours to characterize an 
object region, instead of a conventional circular aperture, 
because we did not want to introduce any possible bias in 
the shape of an object. To define constant metric scale to 
all clusters, we adjusted an extracting threshold in such 
a way that the square root of detected object area times 
a constant was 0.5 Mpc, i.e. const y/area — 0.5 Mpc. We 
chose to use the const =1.5, because the isophotal limit of 
a detected object was best represented by this value. 



3. Morphological measures 

3.1. Centroid & second moments 

Centroid and centered-second moments are computed 
using the first and second order moments of the profile: 
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3.2. Ellipticity 

Ellipticity is simply defined by the ratio of semi-major (A) 
and semi-minor axis (B) lengths as: 



Elli 



l-B/A 



(3) 



where A and B are defined by the maximum and minimum 
spatial rms of the object profile along any direction and 
computed by the formula: 



A 2 = x 2 + y 2 



\ 



B 2 = x 2 +y 2 



3.3. Off-center 



\ 



(4) 
(5) 



The degree of off-center is determined by the distance be- 
tween the centroid and maximum intensity peak: 



Offcen : 



V(* p -x)2 + (y p -y)2 



(6) 



3(A + B) 

where, flux peaks in a pixel at x p , y p . 
3.4. Concentration 

The degree of concentration of the surface brightness profile 
is measured using a method described in Hashimoto et al. 
(1998), and is defined by the ratio between central 30% and 
whole 100% elliptical apertures as: 
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where, is a position of a pixel i in a parameter which 
scales the ellipse, in unit of A (or B), and computed using 
the position angle of each pixel (6i): 
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(8) 



3.5. Asymmetry 



Eh 

ies 



To measure the degree of asymmetry of the profile around 
xtfVfy centroid, an asymmetry index is computed as: 



where Xi and yi are the x-coordinate and y-coordinate of a 
pixel i of value U inside area S of an object. 



Asym = 
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(9) 
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where, 



x2 ieS 



ies 



2/2 



ies 
ies 



(10) 



After testing various centroids, we have chosen to use the 
second order centroid x2 and y2 to make the asymmetry 
measure less sensitive to the very faint outer structure than 
the case using simple x and y. 




4. Uncertainty and systematics 

4.1. Uncertainty 

We applied a Monte Carlo simulation to estimate the un- 
certainties in our measures caused by point sources and 
Poisson noise. For each cluster image, starting from the im- 
age used for real analysis, we added random artificial point 
sources consistent with Chandra PSF and numbers consis- 
tent with the logN-logS given by Campana et al. (2001). 
We chose to use the real image instead of j3 model, since 
many of the clusters were not well descried by the (3 model. 
Poisson noise was then added to the images. We then ex- 
cised the bright point sources again exactly the same way 
as the real analysis, followed by smoothing. For each cluster 
we performed 100 such realizations, and the morphological 
measures were computed for each realization. We then sim- 
ply defined our 1 sigma error for each measure to be rms 
of the distribution of each measure. 

4.2. Systematics 

4.2.1. Exposure time effect 

To investigate the systematic effect of various exposure 
times on the morphological measures, one of the standard 
approaches is to simulate lower signal-to-noise data caused 
by a shorter integration time by scaling the real data by the 
exposure time, and adding Poisson noise taking each pixel 
value as the mean for a Poisson distribution. However, this 
simple rescaling and adding noise process will produce an 
excessive amount of Poisson noise, because of the intrinsic 
noise already present in the initial real data. Meanwhile, 
using a model image with no intrinsic noise, instead of the 
real data, will not have this problem, however, here we need 
to approximate the various characteristics of a model to 
complicated characteristics of a real cluster and this is an 
almost impossible task, particularly for a dynamically un- 
settled distorted cluster. To circumvent this problem, we 
decided to use the real cluster data and employed a series of 
'adaptive scalings' accompanied by a noise adding process. 
Namely, to simulate data with integration t=tl, an original 
unsmoothed image (including the background) taken with 
original integration time tO was at first rescaled by a factor 
Ro/(l-Ro), instead of simple R , where Ro=tl/t0, t0>tl. 
That is, an intermediate scaled image Ii was created from 
the original unsmoothed image Iq by: 



h = h 



Ro 



(1 - Ro) 



(11) 



Poisson noise was then added to this rescaled image by 
taking each pixel value as the mean for a Poisson distribu- 
tion and then randomly selecting a new pixel value from 






Fig. 2. Simulated exposure time (in ksec) vs various mor- 
phological measures for several typical clusters. 



that distribution. This image was then rescaled again by a 
factor (1-Ro) to produce an image whose signal is scaled 
by Ro relative to the original image, but its noise is ap- 
proximately scaled by \/R~o, assuming the intrinsic noise 
initially present in the real data is Poissonian. (The deriva- 
tion of this scaling is described in the Appendix.) Finally, 
the image was smoothed with Gaussian of a— 5" . 

Figure 2 shows the effect of exposure time on each mor- 
phological measure. For each cluster data, we simulated ob- 
servations with several shorter exposure times using the 
method described above, and re-measured our morpholog- 
ical measures for each simulated observation. In Figure 2, 
we plotted the simulated exposure time (in ksec) against 
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various morphological measures of our sample clusters. For 
brevity, we only plotted a handful of typical clusters and 
the data points of the same cluster were connected with a 
line to illustrate the trend. The figure shows that, for all 
of our measures, the systematics caused by exposure time 
are small, demonstrating the relatively robust nature of our 
measures against the exposure differences. The morphologi- 
cal measures are generally constant over a range of exposure 
time, except for some clusters at a very short exposure end, 
where the noise becomes dominant and measures become 
uncertain. 



4.2.2. Redshift effect 

To investigate the systematic redshift effect on the morpho- 
logical measures, we simulated an observation of a cluster 
at a higher redshift than its actual redshift using the real 
data, including the effect of waveband shift, a smaller an- 
gular size of the object (also equivalent of having a bigger 
pixel scale and bigger smoothing scale), and dimming of the 
object signal with respect to the sky. Namely, to simulate 
an observation at new redshift z=zl, at first, we created 
an image of restframe 0.7-8 keV band at original redshift 
zO of each cluster (zl>z0). The image was then corrected 
for detector response and telescope vignetting. Because the 
dimming of surface brightness due to the redshift only oc- 
curs to the cluster signal, and not to the background, an 
object-only frame Io was created from this restframe image 
by subtracting a constant background B. However, to ap- 
proximate the dimming of cluster by the redshift, we can- 
not simply scale Io by 1/(1 + z) because the noise will 
not be correct (we underestimate it). To properly scale Io 
with the proper amount of noise, we employed an adaptive 
scaling technique similar to the exposure time case in sec 
4.2.1. Unlike the exposure time case, however, the intrinsic 
noise contained in Io is not proportional only to Io (it is 
proportional to Io+B, instead, even if the background B is 
already subtracted from the signal). This makes the adap- 
tive scaling more complicated, and we need a pixel-to-pixel 
scaling (or manipulation) rather than just a simple whole- 
image scaling. Namely, an intermediate scaled image Ii was 
created from I by a pixel-to-pixel manipulation: 



h{x,y) 

where, 
Ri 



I {x,y) 2 Rl 



[Jo {x, y) Ri + B - R*( J {x, y) + B)] 



[(l + z0)/(l + zl)] 4 



(12) 



(13) 



Similarly to the exposure time effect in sec 4.2.1, Poisson 
noise was then added to Ii . A new dimmed image I2 whose 
cluster signal was scaled by Ri with respect to the original 
restframe image with proper amount of Poisson noise was 
then created from this noise-added image I ± by a reverse 
pixel-to-pixel manipulation: 



h{x,y) = h(x,y) 



I (x,y)R 1 +B-R 2 [I (x,y)+B] 



h{x,y)Ri 

Finally, adding back the background B gives, 
hix,y) = h(x,y) + B 



(14) 



(15) 



redshift difference between zO and zl. However, this sim- 
ple rebinning again will not correctly reproduce the proper 
amount of noise caused by the angular-size change due to 
the redshift effect. To properly adjust, again underestimated 
noise due to the simple rebinning, the rcbinncd image was 
rescaled by a factor 1/(R|-1), then Poisson noise was added 
by taking each pixel value as the mean for a Poisson dis- 
tribution. The final image was created by rescaling back 
this noise-added image by a factor (R^-l)/^. The factor 
R2 in the denominator is necessary to rcbin the image in 
such a way to conserve the surface brightness. (The deriva- 
tion of these scalings, or manipulations, are described in 
the Appendix.) 





0.0 b L 

0.0 






This dimmed image I 2 should be then rebinncd by a Fig. 3. Simulated redshift vs various morphological mea- 
factor R2 to account for the angular-size change due to the sures. 
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Figure 3 shows the effect of redshift on each morpho- 
logical measure. For one cluster data, we simulated obser- 
vations with several higher redshifts than its original red- 
shift using the method described above, and re-measurcd 
our morphological measures for each simulated observation. 
In Figure 3, we plotted the simulated redshift (s_redshift) 
against various morphological measures of our sample clus- 
ters. Again, for brevity, we only plotted a handful of typical 
clusters and the data points of the same cluster were con- 
nected with a line to illustrate the trend. The figure shows 
that, for most of our measures, the systematics caused by 
redshift are very small. For the concentration index and 
asymmetry index, there may be a slight trend that high 
concentration or high asymmetry objects tend to slightly 
decrease their values as we go to the higher redshifts, thus 
reducing the contrast between distorted and non-distorted 
morphologies. 

4.2.3. Combining the exposure and redshift effects 

Although the exposure time and redshift effects can be 
treated separately as described in sec 4.2.1. and 4.2.2., these 
two effects are often coupled, because low redshift clusters 
are usually observed with shorter exposures than high red- 
shift clusters. Under these conditions, it is often much more 
useful to treat the two effects together, because one can sim- 
ulate even a longer observation than its original exposure 
time, by intentionally reducing the amount of noise to be 
added for the redshift- effect part. With this treatment, to 
compare clusters of various exposure times, we can simu- 
late an observation with 'increased' exposure time for the 
low-z clusters, in stead of standard way of simulating an 
observation with 'decreased' exposure time for the high-z 
clusters, thus we can compare observations of various clus- 
ters without greatly reducing precious signal-to-noise ratio 
of the high-z cluster data. In detail, we modified the final 
steps described in sec 4.2.2. by introducing one more scaling 
parameter R 3 = t2/t0, t2>t0, where t2 is an increased ex- 
posure time, and tO is an original integration time. After the 
rcbinning (i.e. the rebinning after Eq. 15), instead of simply 
rescaling by 1/(R|-1) described in sec 4.2.2., we rescaled the 
image by a factor R 3 /(R2-R3), where (Rj -R3) > 0, namely, 



h 



i?3 



m - #3) 



where, I3 is the intermediate scaled image, and I 2 is 
the dimmed, background re-added, and rebinned image. 
Poisson noise was then added. The noise added image was 
then rescaled back by a factor (R^-Rs^R^ to produce the 
final image whose signal is scaled by R3 relative to I 2 , with 
a proper amount of Poisson noise. (Again, the derivation 
of this scaling is described in the Appendix.) The maxi- 
mum length of integration time we can 'increase' (t2 TOax ) 
is naturally limited by the original exposure time and how 
much we increase the redshift for the redshift-effect part, 
and determined by a relationship: 



R2 — R3 — 0; 



(17) 



which is equivalent to the case when no Poisson noise is 
added after the rebinning described in sec 4.2.2. Thus, 



J L Jl 



limn w 
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Fig. 4. Distribution of the maximum integration time 
(t2 

max) for our sample. 



This t2 max can be also used as a rough estimate of effec- 
tive image depth. The t2 max provides an estimate of the 
image depth much better than conventional simple expo- 
sure time, because t2 max is related to a quantity which is 
affected both by exposure time and redshift, and thus en- 
abling us to quantitatively compare exposure times of ob- 
servations involving targets at different redshifts (e.g. 100 
ksec at z=0.1 and 100 ksec at z=0.9). The distribution of 
t2 max for our sample is plotted in Fig. 4 for the case zl = 
0.9. Several clusters whose t2 max is much bigger than 600 
ksec are not shown in Fig. 4, for brevity. The figure shows 
what the effective exposure times would be, if all clusters 
were at z=0.9. 

Judging from figures 2, 3, and 4 altogether, we have 
decided to modify all of the observations to be equivalent 
of z=0.9 and t=t2 max , to make sure to eliminate even the 
small systematics in Fig. 3, but otherwise to maximize the 
image quality. Meanwhile, four clusters whose original red- 
shifts above z=0.9 were modified only in the exposure time. 
After this stage, to ensure that t2 max is well above the 
low signal-to-noise end, we discarded clusters whose total 
counts are below 300. The resulting sample size after this 
final data preparation is 101. 



5. RESULTS 

5.1. Comparison with other cluster characteristics 



(16) 5.1.1. X-ray luminosity and temperature 



t2 r . 



tm%. 



(18) 



Fig. 5 and 6 show the relation between our morphological 
measures and X-ray bolometric luminosity (Fig. 5), or X- 
ray temperature (Fig. 6) taken from the literature. In Fig. 
5 and 6, we see no obvious trend, demonstrating that there 
is no systematics on our measures due to cluster luminosity, 
temperature, or possibly cluster mass. The lack of trend can 
also mean that more massive clusters do not show more dis- 
tortions, inconsistent with a simple hierarchical structure 
scenario which predicts more massive clusters are younger 
and thus showing more merger activity. Conversely, the lack 
of trend can also mean that interaction of clusters may not 
simply enhance or reduce the cluster global X-ray luminos- 
ity or temperature, although we cannot rule out the pos- 
sibility of such change for a very brief period of time. The 
lack of correlation of X-ray morphology with X-ray lumi- 
nosity or temperature is also reported by Buote & Tsai 
(1996) based on their power ratio analysis. The value of 
the rank-order correlation coefficient, Spearman p, where p 
= f or -f means a perfect linear correlation of a positive or 
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Fig. 5. Comparisons between X-ray cluster morphology Fig. 6. Comparisons between X-ray cluster morphology 
and X-ray bolometric luminosity. and cluster X-ray temperature. 



negative slope, respectively, while p = indicates that two 
variables are uncorrelated, is 0.04, -0.28, -0.02, and -0.06 
for Cone, Asym, Elli, and Offccn, respectively for Fig. 5, 
and -0.35, -0.01, 0.14, and 0.26 for Cone, Asym, Elli, and 
Offccn, respectively for Fig. 6. 

5.1.2. Visual X-ray Classification 

We attempted to compare our measures with X-ray classifi- 
cation by Jones & Forman (JF:1999) in which they visually 
classified clusters into classes: single(S), elliptical(E), pri- 
mary with small secondary(P), doublc(D), offset center(O), 
and complex(C), base on inspection of Einstein X-ray im- 
ages. Unfortunately, JF sample mostly consists of low-z (z 



< 0.05) clusters, and only 15 JF clusters were in our sam- 
ple. Moreover, they were mostly of S or E type, and there- 
fore, a statistically significant comparison was not possible. 
However, a simple comparison between JF S/E type and 
our ellipticity parameter (Fig. 7) already indicates that 'E' 
clusters in the JF sample tend to to show the higher ellip- 
ticity parameters than ! S' clusters, illustrating the consis- 
tency between visual classification and our measures, and 
between Einstein and Chandra observations. 

5.1.3. Beta model profile fitting 

In Fig. 8 and 9, we compare our morphological measures 
with the isothermal /3-model (Cavaliere & Fusco-Femiano 
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Fig. 7. Comparison between our objective "Ellipticity" 
measure and visual X-ray classification (e.g. E=clliptical, 
S=single, P=primary with small secondary) by Jones & 
Forman (1999). 



1976). The single /3-model fitting function is written as 



S(r) = SO 



1+1- 



-3/3+1/2 



B 



(19) 



where SO, r c , (3 and B are the central surface brightness, core 
radius, the outer slope, and a constant background. Fig. 8 
shows comparisons between r c and the measures, while Fig. 
9 shows comparisons between (3 and the measures. 

Again, we only used the values from the literature. For 
our sample, 72 clusters have published values for the single 
beta-model fitting, although some clusters have very dis- 
torted X-ray morphology and cannot simply be described 
by the single or any beta-model. The value of Spearman 
p is 0.79, 0.33, 0.37, and 0.51 for Cone, Asym, Elli, and 
Offcen, respectively for r c plots, and 0.43, -0.01, 0.06 and 
0.17 for Cone, Asym, Elli, and Offcen, respectively for (3 
plots. There is a correlation between Cone (or 1/Conc) and 
r c , while for other measures, there are no obvious correla- 
tions between our measures and (3 model parameters. The 
correlation between Cone and r c shows that clusters with 
small core radii also show high concentrations. Fig 8 illus- 
trates that our robust Cone measure is qualitatively similar 
to the classical morphological analysis based on (3 model fit- 
ting, and that Cone may be used as a robust measure of 
'photon expensive' r c , providing us with a possible alter- 
native to extend the classical radial profile analysis to the 
faint high rcdshift universe. 



5.1.4. Power ratio 

In Fig. 10, we compare our morphological measures with 
the power ratios (Buote & Tsai 1995) taken from the lit- 
erature. The power ratio method calculates the multipole 
moments of the X-ray surface brightness in a circular aper- 
ture centered on the cluster centroid. The powers are then 
normalized by Po- We refer to Buote & Tsai 1995 for more 
detail. In Jeltema et al. (2005), the power ratios are mea- 
sured for 30 out of 101 clusters of our sample. In Fig. 10, 
the left panel shows a comparison between P2/P0 vs. our 
ellipticity, while the right panel shows P3/P0 vs. our asym- 
metry. Error bars for the power ratios are from Jeltema et 
al. (2005) estimated by Monte Carlo simulation for the 90% 
confidence intervals. (The errors from the normalization of 
the background are not included.) Error bars for our mor- 
phological measures, which are estimated by similar Monte 
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Fig. 8. Our X-ray morphological measures compared with 
core radius (r c , in kpc). 



Carlo simulation except that we additionally include the ef- 
fects of point sources, are plotted for the comparison. The 
90% confidence intervals are multiplied by a factor of 1/1.6 
to be roughly comparable to our one sigma errors. The 
P2/P0 vs. Elli plot shows a tight correlation (Spearman p 
is 0.85), which is expected by their similar definitions. Note 
also that the sizes of errors are also comparable. The P3 /Po 
vs. Asym plot shows a weaker correlation (Spearman p is 
0.79). The weaker correlation is probably due to very large 
error bars for P3/P0, (unfortunately much larger compared 
to our robust Asym,) particularly for the clusters showing 
intrinsically high P3/P0 values. Unfortunately, almost all 
of these high P3/P0 clusters are at high redshift (z > 0.4). 
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Fig. 10. Our X-ray morphological measures compared 
with Power ratios. Error bars for the power ratios are from 
Jeltema et al. (2005) estimated by Monte Carlo simulation 
for the 90% confidence intervals. The 90% confidence inter- 
vals are then multiplied by a factor of 1/1.6 to be roughly 
comparable to our 1-sigma errors. 
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are scattered and occupying various places in these morpho- 
logical planes, showing various morphological characteris- 
tics. However, there is a weak to strong trend between each 
set of two measures; the value of Spearman p is -0.62, -0.78, 
0.39 and 0.79 for Conc-Asym, Offccn-Conc, Elli-Offcen, and 
Elli-Asym, respectively. The correlation is relatively strong 
in the Asym-Conc and Offcen-Conc plots. This indicates 
that low concentration clusters generally show high degree 
of asymmetry or skewness, illustrating the fact that there 
are not many highly-extended smooth symmetric clusters. 
Similarly, a correlation between Asym and Elli may im- 
ply that there are not many highly-elongated but otherwise 
smooth symmetric clusters. The correlation between Asym 
and Elli is consistent with the power ratio analysis by Buote 
& Tsai (1996). 



0.6 

beta 



Fig. 9. Comparisons between our X-ray morphological 
measures and /3. 



5.2. Distributions of the measures 

In Table 2, our morphological measures with 1-sigma error 
for the entire cluster sample are listed, while Fig. 11 shows 
the distributions of these measures for the entire sample. 
An average value for each measure is 0.13, 0.16, 0.22, and 
0.05, while median is 0.13, 0.13, 0.21, and 0.04, respectively 
for Cone, Asym, Elli, and Offcen. Figure 12 shows distri- 
butions of cluster morphology in Conc-Asym plane. Each 
cluster point is represented by an X-ray image of the clus- 
ter. The X-ray image is identical to the image used for the 
measurements of morphology. The same Conc-Asym plot 
with 1-sigma error is plotted in Fig. 13, together with other 
measure-measure planes. In Fig. 13, we can see that clusters 



5.3. L-T relation 

In Fig 14, we plot a distribution of cluster morphology in 
the bolomctric luminosity (Lbol) and X-ray temperature 
(Tx) plane using a subset of our sample with available 
literature data. The straight line is a fit from Wu et al. 
(1999) based on 256 low redshift clusters in the form of 
log(Lbol)=2.721og(Tx)-0.92 showing that our sample fol- 
lows the standard L-T relationship (Spearman p — 0.73) 
with some scatter. To be consistent with the fitted line, 
here we used H o =50, fi m =l, Oa=0, cosmology. When we 
subdivide the sample according to their apparent distor- 
tions, using quantitative definition: "distorted" to be Asym 
> 0.12 or Cone < 0.11, and "non-distorted" to be Asym < 
0.12 and Cone > 0.11, and separately plot using open rect- 
angles for the "distorted" clusters and solid ovals for the 
"non-distorted" clusters (Fig. 15), we see that the distorted 
and non-distorted clusters occupy well-defined loci in the 
L-T plane. If we plot distributions of the shortest distance 
from the point to the L-T line (LTD) in the log(Lbol) and 
log(Tx) plane for the two subsets (Fig. 16), with positive 
sign meaning a lower log(Lbol) than the L-T line, we see 
the same trend. For Fig. 16, a K-S test shows the proba- 
bility that the two distributions are drawn from the same 
parent distribution is only 9.56 x 10~ 5 . Meanwhile, we do 
not detect any significant difference in the width of the dis- 
tributions between two subsets (<r~ 0.078 and 0.081, for the 
"distorted" and "non-distorted" clusters, respectively). 
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Fig. 12. Distributions of cluster X-ray morphology in the 
Conc-Asym plane. Each cluster point is represented by a 
X-ray image of the cluster. The X-ray image is identical to 
the image used for the measurements of morphology. 
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Fig. 11. Distributions of X-ray morphological measures for 
the entire cluster sample. 



5.4. Evolutionary effects 

To investigate evolutionary effects, we subdivide our sam- 
ple into low-z and high-z subsets using a border z=0.5. Fig. 
17 shows distributions of Lbol for the high-z and low-z sam- 
ples, while Fig. 18 shows distributions of our measures for 
high-z and low-z samples. A K-S test shows the probability 
that the two distributions are drawn from the same par- 
ent distribution is 0.01, 7xl0~ 3 , 0.50, and 0.15 for Cone, 
Asym, Elli, and Offcen, respectively. Table 3 summarizes 
K-S statistics for various redshift subsets, including the red- 
shift border other than z = 0.5. In Table 3, overall, we see 
no strong evolution in any of our measures for any red- 



shift subset. There is a hint of weak evolution for Cone and 
Asym, in such a way that high-z clusters show more dis- 
torted morphology, consistent with the sophisticated power 
ratio method of Jeltema et al. (2005), but it is unfortunately 
not statistically very significant. Note also that the possible 
weak trend in Cone becomes increasingly insignificant with 
the high-z (z>0.3) vs. low-z (z<0.3) comparison in Table 3. 



6. Summary and discussions 

Using a sample of 101 clusters at redshift z~0.05-l taken 
from the Chandra archive, we quantitatively investigated 
the relationships between the cluster X-ray morphology and 
various other cluster characteristics. The X-ray morphology 
is characterized by a series of objectively-measured simple 
statistics of X-ray surface brightness, which are designed 
to be robust against variations of image quality caused by 
various exposure times and various cluster redshifts. 

We found: (1) Our measures are robust against varia- 
tions of image quality effects introduced by exposure time 
difference, and various cluster redshifts. (2) The distorted 
and non-distorted clusters occupy well-defined loci in the L- 
T plane. (3) Ellipticity and Offcenter show no evolutionary 
effects between high and low redshift cluster subsets, while 
there may be a hint of weak evolutions for Cone and Asym, 
in such a way that high-z clusters show more distorted mor- 
phology. (4) No correlation is found between X-ray mor- 
phology and X-ray luminosity, or X-ray morphology and 
X-ray temperature of clusters, implying that interaction of 
clusters may not enhance or decrease the luminosity or tem- 
perature of clusters for extended period of time. (5) Clusters 
are scattered and occupy various places in the plane com- 
posed of two X-ray morphological measures, showing a 
wide variety of characteristics. (6) Relatively strong cor- 
relations in the Asym-Conc and Offcen-Conc plots indicate 
that low concentration clusters generally show high degree 
of asymmetry or skewness, illustrating the fact that there 
are not many highly-extended smooth symmetric clusters. 
Similarly, a correlation between Asym and Elli may im- 
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Fig. 13. Distributions of clusters in various X-ray measure- 
measure planes with 1-sigma error. 



ply that there are not many highly-elongated but otherwise 
smooth symmetric clusters. 

During mergers, clusters are expected to follow a com- 
plex track in the L-T plane as shown in several numerical 
simulations (e.g. Ricker & Sarazin 2001). Meanwhile, sim- 
ulations by Randall et al. (2002) and Rowley et al. (2004) 
find that even though the temperature and luminosity of a 
cluster varies significantly during a merger, it still follows 
an L ~ T 2 relation. Apart from simulations, the actual ob- 
servational evidence also shows mixed results. For example, 
a study of a 'major merger' cluster A2319 (e.g. O'Hara et 
al. 2004) measuring L & Tx shows little deviation from the 
self-similar L-T relation, which is consistent with the fact 
that bulk properties of clusters either do not change much 
as a result of mergers or change in a correlated way that 



I 

Fig. 14. Distribution of cluster X-ray morphology in the 
L-T plane. The solid line is the fitted line by Wu et al. 
(1999) based on the local cluster sample. 
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Fig. 15. The L-T relation with morphological subsets. The 
solid ovals are the "non-distorted" clusters and open rect- 
angles are the "distorted" clusters. 
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Fig. 16. Distributions of the shortest distance to the L-T 
line (LTD) for "non-distorted" and "distorted" clusters. 



maintains the small scatter of scaling relations. Meanwhile, 
a study of 'merging double' cluster RX J1053. 7+5735 (e.g. 
Hashimoto et al. 2002) shows significant deviation from the 
L-T relation. 

One of the major problem of these studies is their am- 
biguous definition of 'merger', and resulting heterogeneous 
sample of dynamically unrelaxed clusters. These results are 
further complicated by possible bias in measuring L or Tx of 
clusters with various degree of morphological distortion. For 
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Fig. 17. Distributions of Lbol for high and low redshift 
cluster subsets. 



example, Lbol may be systematically underestimated for 
distorted clusters, because the flux from distorted feature 
arc typically diffuse and extended, and may be excluded 
from a luminosity measurement using a finite-sized aper- 
ture. Tx, on the other hand, can be estimated high without 
the extended outer part of distorted clusters. Because it is 
related to the faint signal from the outskirt, these biases can 
be more prominent in the measurements with small effec- 
tive area satellites, such as ROSAT and ASCA. Similarly, 
the presence of a 'cool core' (e.g. Allen et al. 2001) may bias 
low the cluster global temperature measurements and/or 
bias high the cluster global luminosity. If the cool core tends 
to occur in the dynamically relaxed cluster, the observed 
trend in Fig. 14, 15, and 16 may partially be the result 
of this underestimated global temperature and/or overes- 
timated luminosity of the cool core clusters. Indeed, sev- 
eral studies (e.g. Fabian et al. 1994; Allen & Fabian 1998; 
McCarthy et al. 2004) showed that the clusters with very 
large values of 'mass deposition rate' may preferentially lie 
on the high-luminosity side of the L-T relation. The cool 
core clusters may also partially be attributed to the lack 
of difference in the width of the L-T distributions between 
the distorted and non-distorted clusters (Fig. 16), because 
of possible enlarged scatter of the L-T distribution for the 
non-distorted clusters due to the existence of the cool core 
clusters (e.g. Fabian et al. 1994; Markevitch 1998; Allen 
& Fabian 1998; Reichart et al. 1999), while the distorted 
clusters may intrinsically have the large scatter. 

Unfortunately, most of these studies used the mass de- 
position rate to identify their 'cooling flow' clusters. The 
mass deposition rate is, however, mainly based on the infor- 
mation from the surface brightness profile, roughly equiva- 
lent of measuring the morphology of clusters, similar to our 
concentration index. Thus, previously reported behaviors of 
the 'cooling flow' clusters in the L-T plane may be simply 
reflecting the 'morphology effect' which is qualitatively con- 
sistent with our result, and may be showing nothing more 
than the fact that the clusters with high central surface 
brightness behave differently from those with low central 
surface brightness. Indeed McCarthy et al. (2004) reported 
if they restricted themselves only to the presence of central 
positive temperature gradient to select their 'cooling flow' 
clusters, the correlation between the presence of the 'cooling 
flow' and the L-T distribution became unclear. Meanwhile 
studies using simple X-ray core radius of clusters, rather 
than the mass deposition rate, (e.g. Ota and Mitsuda 2004) 
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Fig. 18. Distributions of various morphological measures 
for high-z and low-z cluster subsets. 



reported possible dichotomy in the core radius distribution, 
which may be consistent with the fact that high central sur- 
face brightness clusters behave differently from those with 
low central surface brightness. 

Regardless of the influence of the cooling core, a simple 
comparison between our morphological measures and the 
luminosity (or temperature) (Fig. 5 & 6) shows the lumi- 
nosity or temperature in the literature is not correlated to 
the morphological measures, indicating that we have no ob- 
vious bias in the measurements of L, or Tx for clusters with 
different degree of morphological distortions. We have also 
tested systematics for various satellites and various litera- 
tures, and found no apparent systematic bias among them. 
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No matter what causes the shift in the L-T plane, the 
fact that distorted and non-distorted clusters occupy well- 
defined loci in the L-T plane demonstrates that the mea- 
surements of the global luminosity and temperature for 
distorted clusters should be interpreted with caution, or 
alternatively, a rigorous morphological characterization is 
necessary when we use a sample of clusters with heteroge- 
neous morphological characteristics to investigate the L-T 
and other related scaling relations. There is much left to 
learn about the effect of the cluster dynamical state on the 
bulk properties of clusters of galaxies 

Acknowledgements. Wc thank Peter Schueckcr for useful discussions. 
JPH thanks the Alexander v. Humboldt Foundation for its generous 
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manuscript. 

Appendix A: Derivation of Adaptive Scalings 

A.l. Exposure time effect 

Intrinsic noise No, assuming Poisson noise, contained in the 
original unsmoothed image I (with background) is: 



iVo 



In the intermediate scaled image L after the scaling, the in- 
trinsic noise contained originally in I will be scaled linearly, 
and then new Poisson noise will be added in quadrature. 
Thus the total noise Ni in the intermediate scaled image L 
after adding the new noise is, 



N, = Jlolf/I^ + h 



\lnih+h 



(A.2) 



Meanwhile, in the final scaled image I2, the total noise N2 
will be: 



No 



(A.3) 



Now, the signal-to-noise ratio in the intermediate scaled 
image L (after adding the noise) and the final scaled image 
I2 should be the same, thus, 



J1/JV1 = h/N 2 

h/y/rf/Io + h = h/y/h 

Solving for I\ gives: 



(A.4) 



h 



t2 t2 

1 o I 2 



hh 

(Io - h) 



(A.5) 



However, the scaling from I to I 2 is the scaling factor R 
described in sec. 4.2.1: 



Rq — h/h 
= tl/tO 

thus, L can be obtained from I and R : 

IqIqRq 
1 (I - I Ro) 
Ro 



(A.6) 



Io 



(1 - Ro) 



(A.7) 



A.2. Redshift effect 
A. 2.1. Dimming effect 

Intrinsic noise No contained in the original (this time, back- 
ground subtracted) image Io is not proportional to Io, but 
proportional to I + B, even if the background B is already 
subtracted. Therefore, assuming the noise to be Poissonian, 



No = y/h + B 



(Ai 



In the intermediate scaled image L , the intrinsic noise con- 
tained originally in Io will be scaled linearly, and then new 
Poissonian noise will be added in quadrature. Thus the to- 
tal noise Ni in the intermediate scaled image L after adding 
the new noise is, 



m = yJ(I Q +B)l2/I*+h 



(A.9) 



Meanwhile, in the final scaled image I2, the total noise 
N 2 should include the contribution from the background B, 
even if the background B is remained subtracted, thus 



(A- 1 ) N 2 = y/h + B 



(A.IO) 



Now, the signal-to-noise ratio in the intermediate scaled 
image L (after adding the noise) and the final scaled image 
I2 should be the same, thus, 



/1/JV1 = I2/N2 



h/yJ{I Q +B)l%/l% + h = h/ y/h + B 
Solving for Ii gives, 

h = 



t2t2 
1 1 2 



[P {h + B)-q{Io+B)\ 



(A.ll) 
(A.12) 

(A.13) 



However, the scaling from I to I 2 is the scaling factor Ri 
described in sec. 4.2.2: 



Ri = h/h 

= [(l + zO)/(l + zl)] 4 

thus, Ii can be obtained from I , Ri, and B: 



h 



P R\ 



[I0R1 + B - R\{I + B)] 



(A.14) 



(A.15) 



A. 2. 2. Angular-size effect 

The noise N2 contained in the dimmed, background re- 
added, and rcbinncd image I 2 is proportional to I 2 . 
Therefore, similarly to the appendix A.l., the total noise 
N 3 in the intermediate scaled image I 3 after adding the 
new Poisson noise will be: 



^3 = xli^i/ip + h 



^il/i'2+h 



(A.16) 



Meanwhile, in the final scaled image I4, the total noise N4 
will be: 



N A 



(A.17) 
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The signal-to-noise ratio in the intermediate scaled image 
I3 (after adding the noise) and the final scaled image I4 
should be the same, thus, 



h/N 3 = J4/JV4 

h/y/ri/q + h = h/Vh 

Solving for I3 gives: 
3 (I2 ~h) 



(A.18) 



(A.19) 



Now, in order to scale the final image I4 while conserving 
the surface brightness with respect to the 'prc-rcbinned' 
image (i.e. l' 2 ), 



I± — h I R% 

thus, I3 can be obtained from I 2 and R2: 
r _ h h I R l 

A. 3. Combining the exposure and redshift effects 

Starting from the equation A.19, namely, 



(A.20) 



(A.21) 



h = 



(1'2-h) 



Then, unlike the equation A.20, the final image I4 is now 
related to the original pre-rebinned image I 2 by the new 
factor R 3 : 

R3 = h/I 2 

= TTk = t2/m 



So, 

I4/R3 = I 2 IR2 



(A.22) 



instead of I4 = I 2 /R 2 m the Eq. A.20. Thus, I3 can be 
obtained from I 2 , R 3 , and R2: 



I 2 I 2 R3/R2 
~ h R3/R2) 

{R\ - R 3 ) 



(A.23) 
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Table 2. Morphological measures Table 2 ( continued ) 



Name 



Cone 



Elli 



Asym 



Offccn Name 



Cone 



Elli 



Asym 



Oflfcen 



A3562 


0. 


143±0. 


.010 


0. 


152±0. 


.028 


0. 


159±0. 


.012 


0. 


012±0. 


A85 


0. 


161±0. 


004 


0. 


053±0. 


013 


0. 


124±0. 


005 


0. 


039±0. 


HydraA 


0. 


171±0. 


007 


0. 


134±0. 


.022 


0. 


043±0. 


008 


0. 


003±0. 


A754 


0. 


093±0. 


.004 


0. 


354±0. 


.008 


0. 


323±0. 


.008 


0. 


153±0. 


A2319 


0. 


105±0. 


.006 


0. 


284±0. 


.021 


0. 


250±0. 


.009 


0. 


098±0. 


A3158 


0. 


120±0. 


.005 


0. 


216±0. 


.016 


0. 


075±0. 


.007 


0. 


026±0. 


A3266 


0. 


104±0. 


.005 


0. 


180±0. 


037 


0. 


296±0. 


023 


0. 


106±0. 


A2256 


0. 


098±0. 


007 


0. 


280±0. 


025 


0. 


117±0. 


021 


0. 


129±0. 


A1795 


0. 


138±0. 


003 


0. 


130±0. 


.008 


0. 


187±0. 


004 


0. 


066±0. 


A399 


0. 


100±0. 


.006 


0. 


161±0. 


.026 


0. 


194±0. 


.012 


0. 


102±0. 


A2065 


0. 


123±0. 


.006 


0. 


309±0. 


014 


0. 


207±0. 


010 


0. 


076±0. 


A401 


0. 


118±0. 


.005 


0. 


189±0. 


.023 


0. 


110±0. 


014 


0. 


054±0. 


Zwl215 


0. 


126±0. 


.012 





267±0 


.046 


0. 


110±0. 


.012 


0. 


036±0. 


A2029 


0. 


161±0. 


.007 


0. 


230±0. 


.026 


0. 


070±0. 


.029 


0. 


018±0. 


A2255 


0. 


106±0. 


.010 


0. 


235±0. 


.036 


0. 


142±0. 


025 


0. 


071±0. 


A1651 


0. 


145±0. 


Oil 





205±0. 


.039 


0. 


098±0. 


.040 


0. 


013±0. 


A478 


0. 


153±0. 


.007 


0. 


211±0. 


.026 


0. 


056±0. 


.007 


0. 


009±0. 


RXJ1844 


0. 


158±0. 


044 


0. 


385±0. 


.071 


0. 


202±0. 


.022 


0. 


061±0. 


A2244 


0. 


140±0. 


005 


0. 


094±0. 


014 


0. 


040±0. 


.004 


0. 


008±0. 


RXJ0820 


0. 


122±0. 


.046 


0. 


044±0. 


.074 


0. 


085±0. 


.024 


0. 


017±0. 


A2034 


0. 


124±0. 


.008 


0. 


159±0. 


.021 


0. 


120±0. 


.037 


0. 


028±0. 


A2069 


0. 


127±0. 


013 


0. 


466±0. 


033 


0. 


214±0. 


010 


0. 


036±0. 


RXJ0819 


0. 


151±0. 


.023 


0. 


083±0. 


.073 


0. 


200±0. 


.019 


0. 


050±0. 


A1068 


0. 


161±0. 


.010 





279±0. 


.043 


0. 


076±0. 


.007 


0. 


029±0. 


A2409 


0. 


131±0. 


016 


0. 


129±0. 


.054 


0. 


094±0. 


013 


0. 


041±0. 


A2204 


0. 


180±0. 


.012 


0. 


090±0. 


.043 


0. 


073±0. 


.009 


0. 


022±0. 


Hercu.A 


0. 


147±0. 


.012 


0. 


212±0. 


.065 


0. 


073±0. 


013 


0. 


007±0. 


A750 


0. 


134±0. 


013 


0. 


141±0. 


073 


0. 


074±0. 


.043 


0. 


028±0. 


A2259 


0. 


131±0. 


.019 


0. 


279±0. 


.061 


0. 


126±0. 


.017 


0. 


031±0. 


RXJ1720 


0. 


158±0. 


.006 


0. 


156±0. 


023 


0. 


073±0. 


.006 


0. 


018±0. 


A1201 


0. 


122±0. 


014 


0. 


457±0. 


061 


0. 


292±0. 


011 


0. 


077±0. 


A586 


0. 


144±0. 


.016 


0. 


103±0. 


.056 


0. 


075±0. 


.014 


0. 


018±0. 


A2218 


0. 


117±0. 


.007 


0. 


160±0. 


027 


0. 


103±0. 


.008 


0. 


041±0. 


A1914 


0. 


108±0. 


.006 


0. 


154±0. 


025 


0. 


113±0. 


007 


0. 


089±0. 


A2294 


0. 


132±0. 


.024 


0. 


062±0. 


.074 


0. 


176±0. 


.018 


0. 


064±0. 


A1689 


0. 


161±0. 


.005 


0. 


177±0. 


.019 


0. 


061±0. 


.007 


0. 


020±0. 


A1204 


0. 


156±0. 


011 


0. 


166±0. 


.040 


0. 


049±0. 


010 


0. 


005±0. 


MS0839 


0. 


164±0. 


015 


0. 


064±0. 


.063 


0. 


078±0. 


014 


0. 


007±0. 


A115 


0. 


083±0. 


.006 


0. 


153±0. 


.015 


0. 


424±0. 


.010 


0. 


126±0. 


A520 


0. 


086±0. 


009 


0. 


261±0. 


027 


0. 


153±0. 


010 


0. 


157±0. 


A963 


0. 


151±0. 


008 


0. 


155±0. 


.031 


0. 


080±0. 


008 


0. 


014±0. 


RXJ0439 


0. 


185±0. 


.020 


0. 


113±0. 


060 


0. 


092±0. 


.023 


0. 


018±0. 


A2111 


0. 


119±0. 


024 


0. 


302±0. 


101 


0. 


152±0. 


018 


0. 


101±0. 


A1423 


0. 


168±0. 


.022 


0. 


265±0. 


.086 


0. 


114±0. 


.025 


0. 


017±0. 


Zw0949 


0. 


167±0. 


.014 


0. 


243±0. 


.065 


0. 


062±0. 


.013 


0. 


001±0. 


MS0735 


0. 


206±0. 


.009 


0. 


259±0. 


.025 


0. 


045±0. 


.009 


0. 


001±0. 


A773 


0. 


121±0. 


.010 


0. 


241±0. 


.038 


0. 


102±0. 


Oil 


0. 


030±0. 


A2261 


0. 


159±0. 


.009 


0. 


143±0. 


.038 


0. 


084±0. 


.010 


0. 


016±0. 


A1682 


0. 


103±0. 


.030 


0. 


267±0. 


112 


0. 


349±0. 


.024 


0. 


093±0. 


A1763 


0. 


123±0. 


014 


0. 


302±0. 


050 


0. 


136±0. 


.014 


0. 


036±0. 


A2219 


0. 


119±0. 


.007 


0. 


389±0. 


.014 


0. 


081±0. 


.006 


0. 


048±0. 


A267 


0. 


128±0. 


011 


0. 


266±0. 


041 


0. 


139±0. 


013 


0. 


023±0. 


A2390 


0. 


170±0. 


.005 


0. 


301±0. 


012 


0. 


121±0. 


005 


0. 


031±0. 


RXJ2129 


0. 


173±0. 


.016 


0. 


216±0. 


.064 


0. 


093±0. 


.018 


0. 


025±0. 


RXJ0439 


0. 


147±0. 


.015 


0. 


199±0. 


.047 


0. 


155±0. 


Oil 





037±0. 


A2125 


0. 


140±0. 


021 


0. 


289±0. 


.087 


0. 


247±0. 


014 


0. 


026±0. 


A68 


0. 


131±0. 


023 


0. 
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014 CRSSJ0030 
Oil RXJ1525 

032 MS0451 
023 MS0016 

033 RXJ1121 

021 RXJ0848 
049 MS2053 
Oil RXJ0542 

032 RXJ1221 
025 RXJ1113 

033 RXJ2302 

022 MS1137 
033 RXJ1350 

041 RXJ1716 
029 MS1054 
014 RXJ0152 

019 WGA1226 
014 RXJ0910 

020 RXJ1053 

042 RXJ1252 
069 RXJ0849 

026 

Oil 

036Table 3. K-S probability with various redshift subsets. 
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0.156±0.032 
0.125±0.011 
0.102±0.022 
0.117±0.008 
0.125±0.009 
0.100±0.028 
0.116±0.081 
0.134±0.028 
0.118±0.019 
0.117±0.014 
0.143±0.097 
0.145±0.036 
0.139±0.012 
0.115±0.031 
0.150±0.021 
0.092±0.009 
0.066±0.013 
0.138±0.012 
0.135±0.000 
0.085±0.014 
0.144±0.000 
0.121±0.021 



0.194±0.101 
0.230±0.013 
0.309±0.071 
0.263±0.018 
0.190±0.022 
0.176±0.063 
0.425±0.127 
0.251±0.096 
0.295±0.072 
0.346±0.048 
0.216±0.137 
0.065±0.119 
0.089±0.043 
0.206±0.106 
0.179±0.100 
0.402±0.024 
0.585±0.035 
0.097±0.041 
0.163±0.000 
0.601±0.068 
0.208±0.000 
0.113±0.098 



0.282±0.028 
0.256±0.007 
0.204±0.015 
0.146±0.009 
0.095±0.009 
0.372±0.013 
0.185±0.014 
0.135±0.025 
0.237±0.020 
0.202±0.015 
0.121±0.035 
0.163±0.033 
0.084±0.016 
0.292±0.017 
0.131±0.023 
0.296±0.010 
0.306±0.013 
0.103±0.013 
0.268±0.000 
0.330±0.008 
0.184±0.000 
0.168±0.007 



0.025±0.031 
0.046±0.018 
0.092±0.051 
0.051±0.029 
0.029±0.026 
0.086±0.070 
0.050±0.062 
0.039±0.020 
0.055±0.042 
0.091±0.041 
0.011±0.035 
0.030±0.027 
0.018±0.020 
0.086±0.037 
0.024±0.027 
0.318±0.063 
0.293±0.060 
0.031±0.021 
0.027±0.000 
0.174±0.019 
0.033±0.000 
0.036±0.067 



Q2g Measure 

013 Asym 
019 Elli 
026 Offcen 

028^ 

029 
038 
016 
013 
015 
Oil 
052 
015 
025 
061 
032 
015 
009 
027 
015 
081 
038 
028 
034 
013 
025 
024 
044 
045 
008 
012 
061 
025 
013 
080 
028 
035 
013 
012 
048 
016 
018 
034 
026 
025 
027 
038 
048 
013 



z>0.5 vs z<0.5 z>0.5 vs z<0.3 z>0.3 vs z<0.3 



0.01 
0.007 
0.50 
0.15 



0.03 
0.005 
0.65 
0.14 



0.33 
0.005 
0.50 
0.28 



